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Abstract 

We develop a new renormalization group approach to the large- N limit of matrix models. 
It has been proposed that a procedure, in which a matrix model of size (N — 1) x (N — 1) 
is obtained by integrating out one row and column of an N x N matrix model, can be 
regarded as a renormalization group and that its fixed point reveals critical behavior in 
the large- iV limit. We instead utilize the fuzzy sphere structure based on which we con- 
struct a new map (renormalization group) from N x N matrix model to that of rank 
N — 1. Our renormalization group has great advantage of being a nice analog of the stan- 
dard renormalization group in field theory. It is naturally endowed with the concept of 
high/low energy, and consequently it is in a sense local and admits derivative expansions 
in the space of matrices. In construction we also find that our renormalization in general 
generates multi-trace operators, and that nonplanar diagrams yield a nonlocal operation 
on a matrix, whose action is to transport the matrix to the antipode on the sphere. Fur- 
thermore the noncommutativity of the fuzzy sphere is renormalized in our formalism. 
We then analyze our renormalization group equation, and Gaussian and nontrivial fixed 
points are found. We further clarify how to read off scaling dimensions from our renor- 
malization group equation. Finally the critical exponent of the model of two-dimensional 
gravity based on our formalism is examined. 
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1 Introduction 

There have been a lot of studies on nonperturbative aspects of string theory that have 
revealed that string theory would be nonperturbatively formulated in terms of a matrix 
model or a gauge theory in the large- A limit. In fact, it is well known that one- or 
two-matrix models provide nonperturbative formulations of noncritical strings pQ defined 
in less than one dimension [2]. Furthermore, several large- A gauge theories or matrix 
models are proposed as nonperturbative formulations of critical superstring theories [3]. 
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However, the latter models are incomplete so far because we do not know precisely how 
to take a double scaling limit in which the gauge coupling is tuned to its critical value 
in a correlated way with the large- N limit. In fact, in order to specify the recipe of 
the double scaling limit, we need information on dynamics of theory such as a critical 
exponent. In the case of noncritical strings, we can get it via the orthogonal polynomial 
method for a matrix model, or conformal field theory techniques for continuum (Liouville) 
theory. On the other hand, matrix models or gauge theories proposed as the models of 
critical string theories have more complicated actions with multiple matrices and this fact 
makes these models difficult to solve exactly, by traditional methods like the orthogonal 
polynomials, the Schwinger-Dyson equations, symmetry argument and so on. However, 
we here notice that in order to extract a critical exponent, we do not always have to solve 
a model explicitly. Indeed, we know that universal quantities like a critical exponent can 
be often derived by the renormalization group (RG) approach j3], where we do not have to 
integrate all degrees of freedom in a model. In the case of matrix models, this procedure 
has been materialized in a concrete way which is recapitulated in the next subsection. 



1.1 A brief review of large- N renormalization group 

We review a good example in the case of Hermitian one-matrix model presented by Brezin 
and Zinn- Justin in [5]. For the purpose of examining the critical exponent in the large- N 
limit, they first decompose N x N matrix <pN by 



*>N-1 v 

}N ~\ W t a 



where (pN-i is an (N — 1) x (N — 1) matrix, v is an (N — l)-component vector, and a is 
a real number. Then only v and a, namely one row and one column, are integrated out 
to yield a matrix model of <Pn-i- As illustration, starting from the action 

Sn = Ntr N + 9 -<j>%^ , (1.2) 

then at one- loop level and to the leading order in 1/N expansion, we obtain an action of 
(N — 1) x (N — 1) matrix model, 

SW-i = Ntr N ^ ( Q + £) + 9 -<P%.^j . (1-3) 

The key observation here is that we can repeat this procedure and as such this is quite 
analogous to the usual RG. This motivates us to make a "wave function renormalization" 

^_i=p^_i, p = 1 -^r + °{^)^ ( L4 ) 

so that the action (11. 3ft will have again the standard kinetic term in the (N — 1) x (TV — 1) 
matrix model 

S N ^ = (N- l)tr*_! Q^_ a + j0 v-i) , (1.5) 
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from which we read off the change of the coupling constant as 



9=9 



N 



(9 + 4g 2 



1.6) 



Now it is crucial that since N plays a role of the cutoff in the standard RG, we assume 
the Callan-Symanzik like RG equation 



N 



_d_ 

dN 



F N (g) = r(g), 



;i.7) 



where Fw(g) is the free energy of the N x N matrix model: 



F N (g) = — log Z N (g), 



z N(g) = / d(fi N e 



Sn 



'U 



Namely, ( 11 .7[) prescribes the iV-dependence of F^. By plugging the known exact result, 
we can check that the free energy near the critical point g* = —1/12 in fact satisfies (11. 7p 
and there (3'(g) is related to the critical exponent (string susceptibility) 71 as 



7i 



1.9) 



We expect that at a fixed point of the above RG (11.61) . critical behavior would show up 
and analysis of the RG near the fixed point gives us the critical exponent. This is indeed 
the case. From (II .6p we can read off the /3-function, namely dependence of the coupling 
constant on the cutoff iV as 



-g - V 



1.10) 



and we find the nontrivial fixed point g* = —1/4 and the critical exponent 71 = 2 accord- 
ing to (11.91) . They relatively well approximate the exact values g* = —1/12 and 71 = 5/2. 
Hereafter we call this approach as large- N renormalization group (RG)0 



When we attempt to apply the large-iV RG to matrix models for critical strings, 
we immediately recognize the following issue: interpretation of matrices in these "new" 
matrix models is quite different from that in the "old" matrix models describing noncritical 
strings. Namely, in the former we interpret matrices as carrying information on the 
space-time itself; for example, eigenvalues as the space-time coordinates (of D-branes) 
and off-diagonal components as open strings connecting them. This is in contrast to the 
case of noncritical string where a matrix model is interpreted as a tool for describing a 
discretized worldsheet. Hence according to the spirit of the RG [I], it would be better to 
integrate out off-diagonal components far from diagonal ones, because they correspond 
to highly massive modes in the space-time, not just integrating one row and column as 
in [5]. These considerations tempt us to develop a new large- iV RG in which we assign 
the concept of high/low energy to each matrix element, and we integrate out modes with 

1 More applications of the large- N RG are presented in e.g. [6]. 
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the highest energy. If it is possible, the new large- N RG evidently accords with the spirit 
of the original RG [I] and is expected to have nice properties like locality. This is the 
motivation of the present work. We show that we can in fact define a new large- N RG in 
a well-defined manner based on this idea. The key is utilizing the fuzzy sphere structure 
[7J by which each matrix element carries angular momentum. As a consequence, our 
large- N RG has a nice correspondence with the usual RG in field theory. In formulating 
our large- N RG, we will also clarify in what sense it is "local" and mention its particular 
features. We then analyze fixed points of our RG and propose how to read off a critical 
exponent. 

This paper is organized as follows: in the next section we formulate our large- N RG 
on the fuzzy sphere and derive the RG equation concretely. Fixed points of the RG 
equation are analyzed in section [3j It is desirable that these fixed points can be related 
to critical phenomena which have been observed by numerical simulations [Bl El EE] and 
an analytical study [TT] of fuzzy sphere 4 theory via a matrix model. Our results will 
be compared to these works. Section @] are devoted to conclusions and discussions. In 
appendices we enumerate useful formulas which are necessary in our calculations and show 
some identities satisfied by the fuzzy spherical harmonics. 

2 Formulation of large- N renormalization group on 
fuzzy sphere 

In this section we define a filed theory on fuzzy sphere as a matrix model. Then we 
propose a new large- N RG based on the analogy of the angular momentum realized in 
the space of matrices by the fuzzy sphere. 

2.1 Definition of fuzzy sphere 

We begin with constructing a map from the space of functions on the sphere to that of 
N x N matrices following [7J. A function on S 2 with radius p can be expanded as 



where 6, tp are the polar coordinates, and /, m correspond to the angular momentum, the 
magnetic quantum number (projection of the angular momentum), respectively. YJ m is 
the spherical harmonics 



oo 




(2.1) 



1=0 m=-l 




(2.2) 



which satisfies the orthogonality and completeness conditional 




(2.3) 



2 For later use we change the normalization of Yi m by \/An from the standard one. 
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^En.C^M- '"-^^ . (2.4) 

Yi m (9, ip) can be also expanded as 

Y lm (e^) = p- l J2^ m) n xn --- x ^ ( 2 - 5 ) 

where a;* (i — 1 ~ 3) are the standard flat coordinate of R 3 , and cf™\ is traceless and 
totally symmetric with respect to ix, • ■ ■ ,ii- Y{ m = (— l) m Y/_ m implies that cf™;* = 

In order to define an N x N matrix corresponding to YJ m , we first introduce the spin 
L = (N — l)/2 representation of 577(2), Li (i = 1 ~ 3) which satisfies 

[Li,Lj] = ie ijk L k , (2.6) 

and define 



x l = ocL i: a = p\j N2 — . (2.7) 

Thus N x N matrices x"s satisfy an analogous equation of S 2 : Y^=x^ %2 = P 2 ■ We then 
define aniVxiV matrix TJ m corresponding to Y\ m in (12. 5p by 

t„„ n (2.8) 

h—il 

which we call the fuzzy spherical harmonics hereafter and they are considered to form a 
noncommutative algebra of functions on the fuzzy sphere. It follows immediately from 
the property of cf™: that 

7t = (-l) m T|- m . (2.9) 
Since 7} m is in the same representation as that of Y\ m under the rotation SU(2) of S 2 , 

i 

U(R)T lm U(R)~ l = T lm ,R l mm ,(R), (2.10) 

m'=—l 

where U(R) and R l (R) are the N- dimensional and (21 + l)-dimensional representation of 
R E SU(2), respectively. The N- dimensional representation space of T} m (0 < I < 2L = 
N — 1, — / < m < I) is labeled by an integer s (— L < s < L). From the Wigner-Eckart 
theorem, (I2.10p leads to 

(s\T lm \s') = (-1) L - (^ s l m L s ) R(N,l), (2.11) 
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where R(N,l) is independent of s, s' and depends on the choice of \s). Following the 
convention in [12], we take 

(s\T lm \s'} = (-l) L -° ^) y/(2l + l)N, (2.12) 



so that 



— trjy (TimTi, m ^j — Su'S mm i, (2-13) 
V ^\ T im \s2) (s 3 | T} m |s 4 ) = 6 S1S J S2S3 , (2.14) 



N 

Im 



which can be derived easily from f |2. 1 2[) by using (I A . 2 1) and flA. 3[) . where the trace is taken 
over an N x N matrix. In particular, ( 12.141) leads to 



^J2 tT N(°i T im)^ N (0 2 T/ m ) =tl N (Ox0 2 ), (2.15) 

Im 

- tT N (0!T im 2 T/ m ) = trjvd ti N 2 , (2.16) 



N 

Im 



for arbitrary N x N matrices 0\, 2 . (I2.13P and ( I2.14|) show that 7} m form an orthogonal 
and complete basis of the space of iV x N matrices. Namely, any N x N matrix <ft can be 
decomposed as 



2L I 



= ^^ (plmTlm- (2.17) 
1=0 m=-l 

Comparing (12. ID and (12. 17ft . we find that the space of N x N matrices can be regarded 
as a regularized space of functions on S 2 in which the angular momentum I is cut off 
at 2L — N — 1. It cannot be overemphasized that the space of functions on the fuzzy 
sphere has therefore completely finite dimensions, but that nevertheless they form a closed 
algebra and preserve the rotational symmetry 5*0(3) as in ( I2.10p . 

By taking the large- iV limit with p fixed, it is clear from ( 12 .7\i that 

[x\x j ] = iae ijk x k -> 0, (2.18) 

namely the standard commutative S 2 will be recovered. More precisely, as shown in 
appendix [Bj the structure constant of the fuzzy spherical harmonics 

[Ti imi , Ti 2Tn2 ] = ^ ] fhmi limi fera.3^fam3 ' (2-19) 
Z37713 

tends to that of the usual spherical harmonics in the large- N limit with p fixed. 
From the definition of Ti m ( 12.81) . we find that 

[Li,[Li,T lm ]] = l(l + l)T lm , (2.20) 
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which implies that [Li, [Li, ■]] corresponds to the Laplacian on the unit S 2 . In fact, the 
operator on the unit S 2 L{ = —ie^XjOk satisfies 

C? = -A n , 

Id d Id 2 

An = —A^^-n + -^#2- (2-21) 
sm 9 39 89 snr 9 dip 1 

Hence £j essentially corresponds to the adjoint action of Lj. 

In summary of this subsection, we have constructed the mapping rules: 

1. function — > matrix: 

2L l 2L I 

} ^) = J2J2 <f>i™Yim(8, <p) ->• <t> = &rnT lm . (2.22) 

Z=0 m=—l 1=0 m=—l 



2. integration — > trace: 



^0(^ = 1^0. (2.23) 



Notice that this holds as equality. 
3. Laplacian — > double adjoint action: 

-An<P(9,p) -+ [Li, [Li, ft). (2.24) 



2.2 Field theory on fuzzy sphere 

In the previous subsection we have seen that functions on the fuzzy sphere are defined 
as N x N matrices. Hence considering a field theory on the fuzzy sphere amounts to 
constructing a corresponding matrix model. We have identified iV 2 bases {TJ m } in the 
space of N x N matrices, and it is remarkable that they enjoy the concept of the angular 
momentum (I, m) based on which we will formulate a new RG in the next subsection. 

In the following we restrict our interest to a Hermitian matrix, which corresponds to 
a real function on S 2 . (12. 9p leads to 

4 = <ft 0L = (-ir^_ m , (2.25) 

for (12.171) . As an illustration we consideiEl 

Sn = ^tr N {-^[L^? + + |0 4 ) , (2.26) 

which is a natural regularization of the scalar field theory on S 2 with radius p 

s = / ("2^ {cm ^ + T^ 9 ' ^ + f m ^) • (2 ' 27) 

3 A coupling constant m 2 in front of tr n tfi 2 should not be confused with the magnetic quantum number 
m. Later, the mass m 2 is always accompanied by a subscript TV, and may not be confusing. 
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Note that this action is classically equivalent to ( 12 . 26[) by ( 12 . 23[) and (I2.24p . Expanding <fi 
as in (12.1 7p and using the trace formula (IB. 2[) given in appendix [Bj ( I2.26P can be written 
in the momentum spac^l as 

Sn= J2 2 W + X) + p 2m2 Wlm^lm 
Im 

+^ e n((2/,+i)^) 

(r = l~4) 

x Y(-l)~™ (21 + h U l )( 1 h k \ 1 1 1 U T \\\\ h T k T \ , (2.28) 

y rrii m 4 m J \ — m m 3 m 2 J yL L L j yL L L j 

where the summations effectively run in regions constrained by the triangular conditions 
(mentioned in (lA.lj) ) imposed by the 3j- and 6j- symbols. Hereafter a sum over Z, m 
should be understood likewise and we will not show their range explicitly. 

2.3 Large-iV renormalization group on fuzzy sphere 

In this subsection we demonstrate a new large-iV renormalization group (RG) by taking 
advantage of the fuzzy sphere. As mentioned in the introduction, the main motivation 
to develop the RG on the fuzzy sphere is that there we can introduce the concept of the 
angular momentum in the space of N xN matrices. This contrasts with the original large- 
RG proposed in [5] , where there is no clear notion of momentum. Thus the fuzzy sphere 
enables us to perform a new RG transformation lowering A^: A^ — > N — 1 by integrating 
out matrix elements with the maximal angular momentum. The difference between our 
RG and that in [5] is just the choice of the basis in the space of N x N matrices based 
on which the large-A^ RG is formulated. There would then be no gauge-invariant notion 
of matrix elements of high/low momenta since they are mixed up under U(N) rotations. 
However, introducing the kinetic term breaks U(N) and enables us to define high/low 
momentum with respect to the choice of the kinetic term action. Our formulation then 
manifestly accords with the spirit of the RG [4j. The resulting RG is expected to be 
a local and well-defined transformation and to give rise to a local Lagrangian with nice 
behavior in UV. Notice that for the purpose of the large-A^ RG we have to formulate 
the RG for finite N. The fuzzy sphere which can be formulated with finite A^ fits in our 
purpose, while the noncommutative plane based on the Heisenberg algebra does not, at 
least in a straightforward mannerjf] 

Now we formulate our large-A^ RG. We start from the N x N matrix model as in 
dZ26]) 




(2.29) 



4 Strictly speaking, we are considering an analogy of the angular momentum space, which will be called 
the momentum space for simplicity. 

5 We will mention a serious problem in use of the fuzzy torus in section |4j 
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where we have put the subscript N on the parameters in order to keep track of their RG 
flow. The RG transformation is carried out through the following two procedures: the 
first is coarse-graining, and the second is rescaling. For the first, we integrate out modes 
with the maximal angular momentum, i.e. <pi m with / = 2L and —2L < m < 2L in the 
expansion (12.171) . Then the number of remaining modes will be N 2 — (4L+ 1) = (iV — l) 2 , 
which is exactly the number of components of (N — 1) x (TV — 1) matrix. Thus we would 
reinterpret the resulting action as that of (iV — 1) x (N — 1) matrix. Namely we construct 
a new map from a matrix model with rank iV to that with rank N — 1. Next we make a 
scale transformation p^ — > pn-i in order to recover the original scale, and finally read off 
the change of parameters m 2 N and to rn^^ and g^-\- This is our new RG by utilizing 
the fuzzy sphered In the following, we will present the details of each step in order. 

2.3.1 Coarse-graining 

For notational simplicity, let us divide the space of the angular momentum A as 

A = {(/, m ) | o < I < 2L, -/ < m < 1} C Z 2 , 
Aout = {(l,m)\l = 2L, -2L <m< 2L} c A, 

A in = A \ A out = {(/, m) | < I < 2L - 1, -I < m < 1} , (2.30) 
and define correspondingly 

Thus we formulate the coarse-graining procedure as 

/2L 
II d<PTLe- SN{m "> 9N \ (2.32) 

m=-2L 

where 5 , Ar_i(m 2 v _ 1 , gN-i) is an action of (N — 1) x (N — 1) matrix^ It is important that 
our large- RG formalism respects the rotational symmetry 5*0(3), because it integrates 
out the whole components of one irreducible representation of 5*0(3) which does not mix 
modes with different / as in (I2.10p . 

In order to make calculation of f 1 2 . 3 2 j) tractable, we first note that the kinetic term in 
(I2.29p can be decomposed as 

Q (kin.) _ Pn ( 1 rr n2 1 9 _ o(kin.)in Q (kin.)out 



6 Another RG approach to matrix models based on the configuration space is proposed in |13j . 

7 A similar RG has been proposed in [14] for the purpose of analyzing the noncommutative field theory. 
However, it should be noticed that in our large- N RG we construct an (N — 1) X (N — 1) matrix model, 
not just integrating the out-modes. As a result, we have to take account of renormalization of noncom- 
mutativity as discussed in present subsection, which is in contrast to [14) . where the noncommutativity 
is fixed. 



9 



4 kin - )in = E Iw + ^+pWnW^l, 

(l,m)eA ia 
2L 1 

E 5 (^-l)+M)C^m. ( 2 -33) 



(kin.) out 



m=-2L 

and define 

"1 — r j_ O (kin.)out , / -| — |- , (-.(kin.iout 

II d ^l m e~ S - , (0)o = Y II «™ 0eA > ( 2 -34) 

m=-2L ^ m=-2L 

then our RG equation (I2.32p becomes 

S^-iKr-i^-i) = -logZo + 4 kin - )in -log(e-^° t0 ) o) 

5 (pot.) = ^ triV0 4_ (2 35) 

o / c(p°t-) \ 

Thus the calculation of SN-i(m%_ 1 , gN-i) amounts to evaluating (e~ N ) . For this 



purpose, we classify the potential term by the number of out-modes as 

2 4 
Q ( P ot.) _ Pn9n sr^ T/ 

N ~~ N~ ^ h 

i=0 

2L 

V l = tr N (0 in3 out ) = E A** {r"T 2Lm 



m=-2L 



v 2 = vf + V 2 NP , 



2L 



m,m'=—2L 

2L 



v 2 Np = hr N (r<P mt rr xt ) = \ E ^tL^tL^n {rT 2Lm rT 2Lm> ) , 



m,m'=—2L 

2L 



V 3 = tr N (><r t3 ) = E ^LrnAA^N {rT 2Lmi T 2Lm2 T 2Lm2 ) , 

mi,m2,m3=-2L 
1 2L 4 / 4 \ 

V A = -tr N r nti E II II T ^ > (2.36) 



mi~m4=— 2L i=l \i=l 

where <f) m and out denote a matrix only with in-modes and out-modes, respectively: 

2L 



E <PZT lm , 0° ut = E ^f^m- (2-37) 

(;,m)eA in m=~2L 



10 



The expectation value can be evaluated by using the propagator of the out-modes 

($?m$?m'> = 5 m+m/ (-l) m P N , (2.38) 

where Pn does not depend on m, m! as 

PN = N(N-1) + P %m%- (2 ' 39) 
By using this, we can perturbatively integrate out <p out as 

- log ( e -*r ') =^£<^o-5 2 £ <W5>. + 0(A). (2-40) 



'J 



where the subscript c means taking the connected part. In the rest of the paper, we will 
demonstrate the RG transformation to the second order in g^- Noting the ^-symmetry 
of the original action (I2.29p . we find that 



-,2 



^K-.s.-i) =4 km ° m + ^ (vb + (^)o) 

- 5 {^f) 2 ( W c + <v?) c + (v 3 \ + 2 ^v,). + 2 (v*vi> J 

+ C(^) + const., (2.41) 

where we do not include the expectation values containing only V4, since they do not 
depend on m and then will not affect the renormalization of the parameters m 2 N , g^. We 
postpone carrying out the actual calculations and first explain the rest of the large- N RG 
procedure. 



2.3.2 Rescaling 

The second procedure of the RG transformation is the rescaling. We will see that it also 
involves the renormalization of noncommutativity. After coarse-graining, we should have 
a matrix model of size (N — 1) x (N — 1), say 0, instead of N x N matrix <f) m . To do it, 
we relate m and 4> by 

{-\[L u r}^ = j^j^n-i (-\[L^f) , (2-42) 

which implies that 

C, = iim, (2.43) 

with cf) = m)£A in (fiimTim, where L, is the SU{2) generator (12.61) of spin L — 1/2 = 
(N—2)/2 representation, and T is the fuzzy spherical harmonics with rank N—l. Namely, 
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we fix the wave function renormalization in such a way that the kinetic term is canonically 
normalized. Then it immediately follows that 



On the other hand, it is nontrivial how the tr^r (V n4 j vertex becomes in terms of <p. For 
this purpose it is evident that we have to shift all L's appearing in the 6j-symbols in 
(I2.28P by —1/2 (recall L = (N — l)/2). Here the key recursion relation of the 6j-symbol 
is 115] 



a b n 1 



L L Lj y fl2L + a + l)(2L- a)(2L + 6 + 1)(2L - b) 

-2Ly/(2L + l + l)(2L-l) | L % L li L ii} + V(a+ l)a(6 + 1)6 | 



a 6 / 
LLL-X 

(2.45)" 



We find that the first term indeed provides the necessary factor for the tr/v-i (V 4 ) > while 

the second term does not seem to correspond to natural operation of an (N — 1) x (N — 1) 
matrix. However, as in the usual RG, we are interested in the low energy regime Zj -C L 
(i = 1 ~ 4) in tr Ar (Ti imi T^ m2 Ti 3m3 T Um4 ) and a, 6 in (I2.45j) applied to (I2.28P correspond 
to low energy modes. We then find that the second term is of 0(1/N 2 ) compared to 
the first term, which can be therefore neglected in the large- iV RG discussing corrections 
of 0(1/ N) when N -> N - 1. Substituting ( |2~4"5|) into an explicit form flR3j) . a little 
calculation shows that 

tr* (0 4 ) = (l + 1) tr N ^ (0 4 ) + O (JL^j . (2.46) 

We also rescale p 2 N in our RG transformation. In the ordinary RG on the two- 
dimensional square lattice, the lattice spacing will be multiplied by \[2 after the block 
spin transformation, and we rescale the length scale by l/y/2 in order to recover the 
original lattice spacing. In the case of the RG of field theory, after integrating out high 
energy modes A/6 < p < A (for some 6 > 1) with A the momentum cutoff, we make a 
scale transformation as p — > bp in order to retrieve the original momentum space. Since 
in the case of the fuzzy sphere the maximal momentum is given by the Casimir of the 
mode with 2L, we rescale as 

2L{2L + l) = 2L(2L - 1) ^ 
Pn Pn-i 

namely, in such a way that the maximal momentum will be fixed in a rotational invariant 
manner. Here we define a ratio 

b\ = -fL. (2.48) 
Pn-i 
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Thus b 2 N is given by 

* = i& 2 ' < 2 - 49 > 

by which we make a scale transformation in order to recover the original scale as in the 
usual RG of field theory. Notice that the degrees of freedom in a unit volume would be 
then multiplied by N/ (N — 1), which is consistent with the picture that the fuzzy sphere 
defined by N x N matrices describes N quanta with volume given by a 2 N as in 

°& = wzi&- (2 - 50) 

Correspondingly, it should be noticed that under the scale transformation (I2.47p . a at is 
invariant: 

-l + o(±). (2.51) 



In other words, in the case of fuzzy sphere we take the large-iV limit with the physical 
length scale a 2 fixed, as the continuum limit in the lattice QCD with the physical pion 
mass fixed. Since a is the physical scale characterizing the fuzzy sphere, this limit is 
natural, and accordingly we examine it via the large- N RG which respects a. 

On the other hand, it is known that by taking another large- N limit, we can obtain 
the noncommutative field theory (NCFT) [T7J [18] on the flat two-dimensional plane (see 
e.g. [l2j[T9]). In this case, the large- N limit is taken with the noncommutativity 9 in the 
NCFT fixed as [19] 

2p 2 

N -> oo with 6 = — : fixed. (2.52) 

Hence in order to examine the NCFT via our large- iV RG, we should fix 9^. This implies 
that 

Pn _@n _ On-i _ Pn-i (r> ro\ 

N ~ 2 " 2 _ N-r [Z - b6) 

Namely, in this case we have instead of ( 12.491 ) 

* = j^f- (2-54) 

2.3.3 Change of the parameters 

After evaluating (12.411) and performing the scale transformation, the renormalized action 
can schematically be written as 



c Pn-i , 

ON-x — j-j _ ^ trjv-i 



1 + K N (N, m 2 N , g N ) ) [L h 0] 2 + % ( m 2 N + M N (N, m 2 N , g N ) 



^Pn-i V / z 



i Foi more details, see e.g. [16] 
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+ h ^L N + G N {N,m 2 N ,g N )\tf 

+ (others), (2.55) 

where we have taken account of (12.441) . (12.461) and (12.481) . In general, K(N \m 2 N) g^) ^ 0, 
and a further rescaling of (f> has to be done in order to make the kinetic term canonical. 
The resulting action should be compared with 

1 



S 



N-l 



PN- 



-tr 



7V-1 



[Li, > 



N-l 

Now we read off the change of parameters as 

2 m 2 N + M N (N,m 2 N ,g N ) 



gJV-i | 
4 



(2.56) 



™N-1 



9N-1 



l + K N {N,m 2 N ,g N ) 
9n + G N (N, m 2 N ,g N ) 



(2.57) 



N (l + K N (N,m%,g N )y 

These are RG flow equations for m 2 and g, as N — y N — 1. 

In the last line of (12. 55 p . (others) includes terms of the form which did not appear 
in the original action, induced through the coarse-graining procedure. It contains higher 
order interaction terms tr^-i^ 2 ™ (n > 3) and multi-trace deformations. They would 
induce the RG flow in an enlarged space of coupling constants. However, to the order 
of the perturbation in the present paper, they are not relevant. It also contains various 
"derivative expansions" of kinetic and potential terms, and highly "nonlocal" terms which 
are related to nonplanar diagrams in perturbation theory. In the following subsections, 
we will see how they appear and how to deal with them. 



In following subsections we will concretely evaluate the terms in (12.41 j) and determine 
the corrections shown in (I2.55p . There we will find several quite interesting and charac- 
teristic aspects of our large- iV RG, namely, appearance of "antipode" transformed fields 
and multi-trace operators, and "derivative" expansions in the space of matrices. 



2.4 O(gN) and appearance of antipode transformation 

As a demonstration let us calculate contributions of O(giy) in (I2.4ip . From (I2.36P and 
(I2.38p . we obtain 

<^2 P ) = E (^m^fm'W (r 2 T 2Lm T 2Lm ,) 

m,m' 

= P N ^(-l) m tr w (<T 2 T 2Lm T 2L _ m ) . (2.58) 

m 

Here we note the formula given in (IB. 131) in appendix |B] 

2L 

(-±) m T2LmT 2L - m = (2N - l)l N , (2.59) 

m=-2L 
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(b) 



Figure 1: The planar (a) and nonplanar (b) contribution of 0{g^). The bold line repre- 



sents the propagator of the out-mode (12.381) . 



and get 

(V 2 P ) Q = (2N - l)P N tv N ^ 2 , (2.60) 
which therefore contributes to the mass correction. This contribution comes from a planar 



diagram given in (a) in FigJTJ 

On the other hand, we also have a contribution from the nonplanar diagram shown in 



(b) in FigJB 



m,m 

Pn 



(2.61) 



Hence the nonplanar diagram generates a new operator. In order to elaborate on it, we 
apply the formula in (1B.12j) in appendix |B] as 



J2(-l) m tr N (0 in T 2Lm in T 2L „ m ) = N(2N - 1) E { ^ ^ j 



in j. 

L L 2L 1 N 



(Z,m)gA in 



It is then natural to define an operation on N x N matrices as 

Tim 1 ► ^1 TO = ( — 1) 2]m; 

thus 



y n (-i)^ m ) . 

(2.62) 
(2.63) 
(2.64) 



(Z,m)eA in 



(Z,m)eA in 



Although this transformation looks strange in the space of matrices, it has a natural 
interpretation on the fuzzy sphere. In fact, it corresponds to a discrete transformation on 
S 2 



Y lm (9,<f) 



[-l) l Y lm (0,<p)=Y lm (ir-9,<p + ir), 



(2.65) 
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namely, it transports a field to the antipode point on S 2 . (I2.64p is therefore a natural 
counterpart on the fuzzy sphere of this transformation and hereafter we call it antipode 
transformation. Using this, ( I2.6ip becomes 



(Vr) = ^N(2N-l) £ {lL2l}^(^)- 



(2.66) 



(Z,m)eA ; 



Furthermore, we can make a sort of "derivative expansion" of this expression. In fact, 
from (|A.9[) the 6j-symbol in (I2.66P has a large- L (thus large- N) expansion^) 



L L I 
L L2L 



-r. 



2L + 1 
1 

2L + 1 



P, I 1 
1 



1 



L + l 
1(1 + 1) 



2{L + 1) 



O 



V- 



(2.67) 



to give 



Pi 



N 



(2N - 1) 



tr 



A' 



W - 



N 



:tr 



N 



+ o 



1 



(2.68) 



where we recall L = (N — l)/2 and fl 2 . 2 1) . It is quite interesting that owing to the fuzzy 
sphere structure, the space of N x N matrices admits the derivative expansion of an 
operator. Here higher derivative terms come from the /-dependent terms of the expansion 
of the Legendre polynomial in ( I2.67P in terms of 1/L and are suppressed by 1/N as such. 
Notice that this expansion should be in terms of the Laplacian (I2.20p . because the fuzzy 
sphere preserves the rotational symmetry SO '(3), which is respected in our large- iV RG 
formalism. Thus we find that the operator on the right-hand side in ( I2.66P contains 



tr^v ^0 m m "J and its derivative^. We point out here that new operators that are not 

contained in the original Lagrangian ( I2.29P are thus generated in our RG, but that they 
may be regarded as small corrections. For example, the first term in (I2.68P agrees with 
the overlap between a wave function and its antipode 



N 



:tr 



N 



dVL 
Air 



E 

(Z,m)eA in 



-1)' 



Im I 



(2.69) 



which is always smaller than 4?trAr (</> m2 



iV-iv \y j ^ 2~2 { i >m )eM n l^ml 2 in (ESQD coming from the 
planar diagram. In fact, it is easy to check that for the Gaussian wave function (under the 
stereographic projection) the former is exponentially smaller than the latter with respect 



9 Throughout this paper, we concentrate on the case with the "in" fields carrying low enough momenta. 
Therefore in the expression of the perturbation theory, we write errors simply as 0(1/L 2 ) or 0(1/N 2 ), 
though it could also depend on the momenta of associated "in" fields. 

10 At a glance, the antipode transformation seems to correspond to exchange of two matrices. However, 
as shown in appendix [D] it is related to the reverse of the ordering of all matrices inside a trace. 
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to the inverse square of the width. However, since in the RG we are interested in low 
energy modes, it is not obvious that the overlap is small enough to ignore. Thus in the 
fixed point analysis given in section [3] we will concentrate on the flows of m? N and 
up to 0{g 2 N ) and simply assume that to this order, interactions involving (p mA would not 
affect them considerably, namely, would not modify ( I2.57P so muchj^j] Here it is important 
to note that as we see by comparing the coefficients in ( 12.601) and (12.681) . the nonplanar 
contribution is not suppressed by 0(1/ N) in itself compared to the planar one, but the 
momentum dependence (in the present case (— 1)' in the antipode transformation) can 
make it suppressed. Recall that this is also the case with the NCFT, where nonplanar 
diagrams are in general dropped in a high energy regime by the oscillating phase depending 
on momenta. 

2.5 0(g 2 N ) and appearance of double trace operator 
2.5.1 Vertex correction 

Among contributions of 0(g 2 N ) shown in the second line in f )2.4ip . 4 vertex correction 
originates only from (V 2 2 ) c . Thus in this subsection, we focus on this part and present 
only the results of the calculation of the other contributions. The details are given in 
appendix O 

Since V 2 = V 2 P + V 2 NP as in ATM , let us begin with (v 2 p2 ) . From f l2~36l) . 



n> 2 

x tr N T 2Lmi T 2Lm2 J tr N ^0 in T 2Lm >T 2Lri 



mi,m,2,m' 1 ,m' 2 



Pn J2 (-l) mi+m2 (trv {^ 2 T 2Lmi T 2Lm2 )ti N (0 in2 T 2i _ m2 T 2i _ mi ) 

rn 1 tyi n \ 

+ tr^v T 2Lmi T 2Lm2 ^j trjy (4> m T 2L _ mi T 2L _, r , 



mi,m2 



(2.70) 

The first and second term correspond to the planar and nonplanar diagram shown in 



Figj2] [(aj| and (b) , respectively. Important observation here is that double trace terms as 
above are in general generated in the large-iV RG. At first sight, they seem to be new 
interactions arising in our RG and do not play any role in renormalization of parameters 



11 'A 

Here we should notice that it is likely that an interaction with <p ln originating from a nonplanar 
diagram is not irrelevant in itself near a hxed point. In fact, it is shown in [19] that a contribution from 
the 1PI nonplanar two-point function r^j* Janar (corresponding to (V^ p ) with all modes running the 
loop integrated) gives rise to the IR singularity in the large- N NCFT limit (|2.52[) and hence becomes the 
origin of the UV/IR mixing. Since the large- N limit should be captured by a fixed point in the large- N 
RG, this suggests that if an interaction involving (f> ln has something to do with r iio ; nplanar which becomes 
quite large in the IR in the large- N limit, it would be rather relevant around a corresponding fixed point. 
We will further make a few comments on such an interesting aspect in section |4j 
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(a) 



(b) 



Figure 2: The planar (a) and nonplanar (b) contribution in (V^ 2 ) . The bold line 



represents the propagator of the out-mode (12.381) . 



Pn, m 2 N and g^ in the original action (12. 29ft . But this is not the case. Recall that in the 
usual RG of the 4 scalar field theory we also get (in the configuration space) from the 



diagram in (a) in Fig. [2] 



dxdyA(x — y) 2 4>(x) 2 (fi(yy 



(2.71) 



where A(x — y) — (<f)(x)(f)(y)) is the propagator of the scalar field. Namely, bi-local 
interactions are generated in general in the RG corresponding to the double trace term in 
(I2.70p . However, we know that if the propagator A(x — y) is for a highly massive mode of 
<ft, it would rapidly damp as x and y are separated. In such a case the derivative expansion 

as 



dx 



dzA(z) 2 4>(x) 4 + 



dz—A(z) 2 
2 v ; 



<t>(x) 2 d 2 J(x) 2 + 



(2.72) 



is valid at least for a low energy mode of of interest in the RG. Namely, bi-local 
interactions generated by the RG can be decomposed as a sum of local interactions and 
there higher derivative interactions are suppressed if we integrated out sufficiently high 
energy modes. In other words, the RG by integrating only high energy modes only 
generates local interactions. Notice that the first term in (I2.72p contributes to 4 vertex 
correction. 

Hence in our large- N RG we also have to make a "derivative" expansion as in (I2.72p 
in the space of matrices and to read off what kind of "local" interactions arise as a result 
of the RG. Interestingly enough, it is indeed possible in our formalism because our RG is 
based on the fuzzy sphere structure and we in fact integrate out modes with the maximal 
momentum. As illustration let us consider the planar contribution of (v^ 2 ^ given by 
the first term in (12 .70 p . By using the fusion coefficient (IB.4[) with (IB.5|) in appendix |B| 



V, 



P2 N 



planar 



-P 2 (4L + 1) 2 N (- 1 ) 



mi+m 2 



mi ,rri2 



J2 V(2i + i)(2r + i)(- 

(l,m),(l',rn') 



\l+m+l'+m' 
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2L 2L I 
mi ni2 —m 



x tr N m T lm trjv 



2L 2L V 

-7712 —77l\ —77l' 
J- I'm 1 



2L 2L I 
L L L 



2L 2L V 
L L L 



(2.73) 



Then we can take the summation over m\, first by the formula ( 1A.2I) in appendix [A] 
to get 



V, 



F2 



planar 



~-P 2 N (2N - l) 2 iV^ 



Irn 



2L 2L I 
L L L 



trjv 4> in T im tr N 



1 lm 



. (2.74) 



Now we make a "derivative" expansion. In the large- N RG we are interested in sufficiently 
low energy modes with ! C L as in the usual RG in field theory, and then ( 12.131) 
tells us that the internal momentum I in (I2.74p should be also sufficiently smaller than 
L so that the traces will not vanish. Thus as far as <f> m only with low energy modes is 
concerned, we can evaluate the 6j-symbol in (I2.74p under I ^ N by using flA.9j) as 



2L 2L I 
L L L 



-11 



Pi 



v/(4L + l)(2L + l) \ V 2L + 2 



2L + 1 



y/(2N-l)N 

Plugging this into (I2.74p . we obtain 



1(1 + 1) 
AN 



(2.75) 



P 2 



planar 



P 2 N (2N - 1) Wtr* {^T lm ) tr N (f^T^ 

lm 



~ 2N tlN 



[Li,[Li,T lm ])) tr N U ta2 T* 



' lm 



(2.76) 



Then we see that the first term in fact becomes the single trace by virtue of the complete- 
ness of T; m (I2.15p . It is rather technical but interesting that in the second term we can 
also take the summation over I, m by applying (12.151) . because by the "partial integration" 



tr 



N 



[Li, [L^ Ti m ]] ) tr 



1 lm 



tr 



N 



Tim ) t r iV 



L in 2 rTif 

; 1 lm 



(2.77) 



As a consequence we have 



planar 



Pjj{2N - l)N 



tr 



N 



tr /v 

2N N 



Li, 



+ o 



1 

iV 2 



(2.78) 
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Namely, we have the vertex correction of 0{g 2 N ) from the planar diagram. Detailed 
calculations presented in appendix O show that the other contributions in (V^ 2 ) c always 



generate 



vertex corrections with both d> m and 



nonplanar 



Pjj{2N - l)N 



trW0 in V nA2 



2N tTN 



Li. 



42 
1^1 



[v 2 p vr) c 



P%(2N-1)N 



tr 



N 



V Np2 



\p 2 (2N-l)N 



2N 



trjy ( in in V n inA 



Li. 



+ 



1 

N 2 



1 

iV 2 



(2.79) 



where we notice that all derivative corrections come in as the Laplacian which reflects the 
fact that our large- N RG as well as the fuzzy sphere respects the SO (3) symmetry, as 
anticipated below (I2.68p . Therefore the RG again generates new interactions containing 

■ A 

both (f) m and m that are not present in the original action. For the same reason as in 
(I2.68P we regard them as small corrections and negligible in discussion of the RG of the 
parameters p 2 N) m 2 N and g n in the original Lagrangian (I2.29p . 



2.5.2 Mass correction 

Let us examine the other terms than (V 2 2 ) c in the second line in ( I2.41j) . It is obvious that 
they contribute to correction in the quadratic order of m . The results are summarized 
as: 

1. (V 2 ) c vanishes in sufficiently low energy regime where the momentum I of in is 
much smaller than the cutoff 2L. 

2. (V 2 ) c is exponentially suppressed by L when L ^> 1. 

3. (ViV3) c = by the momentum conservation. 

4. (V 2 p V4) c ~ Pfj2(2N — l) 2 Nti N (V n2 j , which arises from the planar diagram. 

5. (y^ p V^ c again provides derivative expansions of terms with both the usual and 
the antipode fields. 
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Therefore only (VfV^ part contributes to the mass correction. The details of the cal- 
culations are shown in appendix IC1 



2.6 Flow equations 



So far we see that when m is in a sufficiently low energy regime and L > 1, ( I2.4ip 
becomes 



Sat-i — tr 



inl2 , PN m N ,in2 , Pivfi'jV / 1 , in 4 



2Ar L t.V J 2]V V 



iV U 



+ P N {2N - 1)0" 



+ 



N 



trjy 



1 



2 L 



" :|2 + ;(< + ^ J Bi(iV,m^)-p^ J B 1 (iV,m^)S 2 (iV,m^ 4 '" 2 



(g N -PN9 2 N B2(N,m 2 N ))^ 



+ 



(2.80) 



where 



Bi(N, m 2 N ) 
B 2 (N,m 2 N ] 



2{2N - 1)P, 
2(2N - l)Pj 



N, 



JV) 



(2.81) 



with Pjy given in (I2.39p . In the above, ■ ■ • denotes the terms including <ft m independent 
constants, 0(g%), antipode fields, and 1/N suppressed derivative expansions. They are 
all omitted in the following RG analysis. Now the RG flow equations ( 12. 5 7ft become 



m 2 N _ x = b 2 N (m 2 N + g N Bi{N, m 2 N ) - p 2 N g 2 N B x (N, m 2 N )B 2 (N, m 2 N )) 
9n-i = b 2 N (g N - p 2 N g 2 N B 2 (N } m 2 N )), 



(2.82) 



with keeping 0(g 2 N ) terms. 

In the next section we analyze the RG transformation (12.821) . Before closing this 
section we comment on a couple of features of our large- N RG equation in order. 

1. The mass corrections of O(gisr), 0(g 2 N ) and the vertex correction of 0(g%) arise 

from the planar diagrams in (V 2 P ) Q , (V^V^, and (v^ 2 ^ , respectively. Nonplanar 

diagrams generate terms with the antipode fields, which are neglected in the current 
analysis. 

2. No kinetic term correction appears. This reflects the fact that the result of (V 2 P ) 
in ( I2.60P is exact and does not have any derivative corrections. 

3. Since our large- iV RG is constructed by integrating out the modes with the highest 
energy, it should be in a sense local as in the case of the usual RG. In fact, all 
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corrections carry the propagator Pn given in (I2.38p . which is highly suppressed for 
large- N. Furthermore, our RG admits the "derivative" expansion as mentioned in 
section |2"31 by which the double trace corrections can be written as a sum of the single 
traces. This fact also reflects locality of our RG, because the reason why we can 
utilize the asymptotic formula of 6j-symbol as in (I2.75P is that we have integrated 
the modes with / = 2L and hence its upper entries become the cutoff 2L itself. 
Accordingly, we have observed that the "derivative" corrections are suppressed by 
1/JV as in (EZHI). 



3 Fixed point analysis 

In this section we examine Gaussian and nontrivial fixed points of our RG transforma- 
tion (I2.82p and discuss their properties. We also determine the scaling dimensions of 
small perturbations around these fixed points from (12.821) . First we discuss fixed points 
corresponding to field theories, and then turn to the two-dimensional gravity originally 
considered in [5|. Results in this section are compared to preceding works [8j[9|[10j[llJ, 
regarding nonperturbative studies of 4 theory on the fuzzy sphere via the matrix model 
by different methods. 



3.1 Gaussian fixed point 

It is evident that m* = = is indeed the fixed point (Gaussian fixed point) of the RG 
transformation (12.821) . Near the Gaussian fixed point, we linearize (12.821) with respect to 
small perturbations m 2 <C 1 and g <C 1 as 

m N-l = b N { m N + 9nBi(N)) , 

9n-i = b 2 N g N . (3.1) 
Here we expanded Bi(N,m 2 N ) and £? 2 (iV, mfy) as 

B l {N,m 2 N ) =B l {N) + B hl (N)m 2 N + • • • , B 2 (N,m 2 N ) = B 2 (N) + B 2il (N)m 2 N + • • • , 

(3.2) 

where 

p (i\t\ _ 2(2./V" — 1) 2(2JV-1) 

Bl{N) ~ N(N-l) ' B2{N) ~ N\N-lf ■ (3 - 3) 

In order to diagonalize the RG transformation near the Gaussian fixed point, we assume 

the eig envector as m 2 N — m 2 ((?iv), where m 2 (^7v) = A^gN + B^g 2 ^ H ■ We then need to 

solvql 2 ] 

m 2 N _ x - m 2 (g N _i) = b 2 N (m 2 N - rh 2 (g N )) , 



12 As the readers can see, the coefficient An in m 2 (g) depends on the step of RG transformation, namely 
varies as N — > N — n. Therefore m 2 (g) is slightly different from the critical line of the RG flow. The 
leading behavior of An shown here approximates the slope of critical line in the vicinity of the fixed 
point. 
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9n-i = b 2 N g N . (3.4) 
These equations determine fh 2 (g^) as 

m 2 (g N ) = - (41ogiV) g N + O f^j . (3.5) 

Thus rh 2 approximately give a critical line of the RG flow, and then we find that the 
critical line is vertical in m 2 — g plane in large- N limit. This feature originates from 
the fact that m 2 N and have the same scaling dimension and thus it is peculiar to two 
dimensions. In fact, this is also the case with the RG of the ordinary two-dimensional 
field theory where the critical line becomes vertical. This vertical line can be derived in 
a different manner. In order to regularize the slope of the vertical line, we introduce an 
extra tiny iV-dependence c(N) ~ N~ £ to g n in f !2.29|) 

c Pn+ ( 1 \t a 12 , m N,2 , c(N)g N A 

S N = —tr N ^-^[L U (j) N } +—<P + ~ A J, (3.6) 

and repeat our large- N RG to obtain the RG equation 

m 2 N ^ = b 2 N (m 2 N + c(N)g N B 1 (N,m 2 N ) - c(N) 2 p 2 N g 2 N B 1 (N,m 2 N )B 2 (N J m 2 N )) , (3.7) 
9n-x = ^X/ n (9n ~ P 2 N c(N)g 2 N B 2 (N, m%)) . (3.8) 

Now critical line ml(g) near the Gaussian fixed point can be determined as 

9/ . 12(2JV- 1) 4 . , 

m 2 M = — £ \_ x } 9^--9 (N^oo). (3.9) 

Thus we also find the vertical critical line. Compared to the case of the ordinary field 
theory, we see that the introduction of c(N) ~ iV _e is an analog of considering (2 + e)- 
dimension. Furthermore, (13. 7p and (13. 8p tell us that the eigenvalues of the RG trans- 
formations are b 2 N for m 2 N and J^_ x \ b% fc> r gN- Now we propose how we read off the 
scaling dimension in our large-iV RG. Recall that we make the scale transformation by 
= Pn/Pn-1 m or der to compensate the maximal momentum integrated out in the 
large-iV RG. Therefore, the scaling dimension should be given by how many powers of 
h n an eigenvalue of the RG transformation has, as in the standard way to extract the 
scaling dimension in the usual RG. In other words, since = 1 + 0(1/N), we can read 
off the scaling dimension by observing deviation of 0(1/ N) from 1 in an eigenvalue of the 
linearized RG transformation^ Applying this method to (13. 7p and (13. 8p . we conclude 
that our Gaussian fixed point is an UV fixed point and the scaling dimensions of m 2 and 
g are both 2 for infinitesimal e. It is in fact consistent with the result from the mean-field, 
or the Landau theory, as it should be near the Gaussian fixed point. 



13 Since we have just changed N by one in the large- N RG, the eigenvalue must be 1 in the leading 
order in the 1/N expansion. 
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The critical line found in present analysis also have been seen in numerical works 
[HI Ellin]- The uniform-disordered phase transitions are reported on this critical line. Our 
results shows that the slope of the critical line diverges in large- N limit, while [H] reports 
weak iV- dependence of the critical line. This is because our slope of the critical line is 
logarithmically divergent, and thus it is hard to see the vertical of the critical line if N is 
not so large. 

It is worth noticing that (I3.7j) and (13. 8 ft also imply the Gaussian fixed point in the 
NCFT, where the scaling dimensions of m 2 N and gpj are again 2, because the difference 
between the field theory on the fuzzy sphere and the NCFT is just in the iV-dependence 
of 6jv itself as in (I2.49j) and (j2.54|) . namely in a quantity which should be fixed in the 
large- N limit. It would be intriguing aspect of our large- N RG formalism that we can 
discuss a field theory on the fuzzy sphere and on the noncommutative plane in a unified 
way. 



3.2 Nontrivial fixed points 

Next we search for nontrivial fixed points. If there is a nontrivial fixed point, we would 
have nontrivial wave function renormalization there. In such a case if we normalize the 
kinetic tern canonically as in ( 12.421) . it would give rise to nontrivial iV-dependence of 
as in (13. 6p . Here we assume general leading iV-dependence of c(N) in (I3.6P as c(N) = cN a 
and fix the value of a in such a way that there exists a fixed point of O(N ). Near such a 
fixed point, if any, the perturbative expansion in terms of can be in danger due to higher 
order interactions. We observe that since p 2 N c{N)g^P^N is the loop expansion parameter, 
p 2 N c(N)g N P^N <C 1 is required for the loop expansion to be a good approximation^- 
With this condition the third term in the left hand side of (13. 7p is subleading. Therefore 
up to O(gff) the RG equations (13.71) and (13. 8p can be recast into 

™N-i = b% (m 2 N + c(N)giyBi(N, m 2 N )) , (3.10) 
1 _ 1 

We find a fixed point of (13.101) and (13.111) analogous to the Wilson-Fisher fixed point 

2 _ N(N- l)^(^y^-l) _ N 2 (N-l) 2 ^r ) b%-l ( b%-l V 
~ & =583^-1 ' * 2(2AT-1) c(N)p 2 N {^ b% _ 1 )' 

(3.12) 

We then confirm that the loop expansion parameter p 2 N c(N)g]siP^N is small around this 
fixed point as a consistency check. Using these fixed point values, we evaluate 

Ao(N )g .PlN = (-^^ - l) . (3,3) 

14 Here we do not take further possible suppression factors due to 6j-symbols into account. 



9N 



+ c(N)p 2 N B 2 (N, 



m 



N, 



(3.11) 
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For arbitrary a, this is of 0(N~ l ) and negligible in the large- N limit. Therefore our fixed 
point ( I3.12j) is stable against higher order loop corrections, and this also suggests that the 
fixed point will not receive much effects in general from higher order interactions induced 
by the RG transformation. Such feature can be understood as the locality of our RG 
transformation itself. 

We linearize (I3.10P and (13. lip around the fixed point. The linear perturbations obey 



hi (2 - ^Tt/n) Sm% + (1 - b%) ^6g N} (3.14) 
1 / 1 \ (1 - b 2 ) m 2 

6 ^-i = ^r^J 9N + 1 " J b 2 N N ( 2iv-iV(iv) ^- (3 - 15) 

c(N-l) U N \ c(N-l) U N/ N \ 

The scaling dimensions are determined from the eigenvalues of these linear perturbation 
equations. Note that off-diagonal terms in these equations are of 0(N~ l ). Thus to the 
leading order, the eigenvalues for 5m 2 and 5g are 

^ ,.2 \ r . ( C(JV) ^ 



with 0(N~ 2 ) corrections. As concrete examples, we discuss scaling dimensions of these 
linear perturbations in both cases of fuzzy sphere and NCFT, separately. 

fuzzy sphere case: Substituting p 2 N ~ N 2 a 2 /4 where a 2 = a 2 N is fixed in the RG, we 
find 

where we have used (I2.49|) . We choose a = to have of 0(1), then (m^g*) — > 
(— 2—5). Since b N ~ 1 + 1/N from (12.491) . the scaling dimensions of the linear pertur- 
bations are obtained from coefficients of 1/N terms of the eigenvalues (13.161) . They are 
for 5m 2 and —2 for 5g. Namely there are one marginal and one irrelevant perturbations 
to this order in this large- N limit. 

NCFT case: p 2 N = N9/2 from (12331 yields 

m ^——2' v-^^FJ^W {N ^ X) - (3 ' 18) 

where we have utilized (I2.54p . We set a = 1 to have of 0(1), then (ml, g*) — > (— -^-). 
Although the position of ml goes to negative infinity as N goes large, the condition for 
higher loop corrections to be suppressed is still satisfied. Note that now the rescaling 
factor is given as b 2 N = 1 + -k. Thus the scaling dimensions of linear perturbations can 
be read off as —2 for Sm 2 and —4 for 5g. This IR fixed point is stable in this large- N limit. 
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large p 2 N limit: When p 2 N m 2 N 3> A^ 2 , the factor from the propagator Pjv in Bi(N,m 2 N ) 
and B 2 (N, m 2 N ) can be approximated by p^m^ 2 . Then we easily observe from (13. lip how 
the position of the fixed point depends on iV as (for simplicity, we choose c(N) = 1) 

N 1 PN 

This suggests a phase transition which occurs on such points. Indeed, the so-called 
disordered-matrix phase transition is observed on such points in numerical simulations 
[EJ [9j [10] . In [UJ this phase transition is explained by a topology change of matrix eigen- 
value distribution. In these previous works, the matrix model action is given as 

4:71 

S = — tT(<t>[Li, [L h 0]] + rR 2 ^ + \R 2 ^). (3.20) 

The relation between (R, r, A) and our (p N ,m%, g^) is determined as piy = R, m 2 N = r, 
9n = 2T- I n [ 1 1 J , noncommutativity is fixed as R 2 = thus we adopt b 2 N = -J^j as in 
( I2.54p for a comparison. Then the equation (13.191) becomes 

^ = ±7! = ±0 - 564 ir <^°°>- < 3 - 21 ) 

On the other hand, the position of the critical point on A = 1 line is calculated numerically 
in [8], and also estimated by an analytical method in [11]. Their results are 

, ,_, r 0.56 
observed m [8] : 



N R ' 

, • rr^n T 3 1, 0.846 

estimated in 11 : — = ± — -=— = ± . 

J 2^R R 

Our result (I3.2ip have a good agreement with the numerical simulation^!. It strongly 
suggests that our formula (I3.19P describes the disordered-matrix phase transition points. 
Analysis of critical behavior of various order parameters will provide further evidence. 
And the phase structures of this matrix model are further investigated numerically in 
[TO] . It would be interesting to understand these structures by our RG methods. 



3.3 Two-dimensional gravity 

As a final application of our large- A" RG, in this section let us consider the two-dimensional 
gravity. For comparison we start from the same action as in the original large- A" RG in [5] 
and apply ours, namely, perform the large- A" RG by a different basis of N x N matrices 
from that in [3]. Then we compare our results with those in [5J. 

We begin with the one-matrix model (II. 2p defining the two-dimensional quantum 
gravity 

Sn = Ntr N + , (3.22) 

15 If we adopt b% = as in (pQ9l) . the number in front of ^ in (pT2T)l differs by 2" 1 / 2 . 
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here dropping the kinetic term | [Li, [Li, fa]] would be essential for the two-dimensional 
gravity, because by doing so we have the U (N) gauge symmetry in (I3.22p which would be 
somehow related to the diffeomorphism invariance. Then it is straightforward to repeat 
our formalism. Two important differences from before are that now we fix the wave 
function renormalization as 

Ntr N (> 2 ) =(N- l)t rjv _ a (?) , (3.23) 
instead of (I2.42p . and that the propagator becomes 

Pn = ^. (3.24) 
With these in mind, it is easy to derive the RG equation for g N again as 

9N-i = 9N ( 1 - |f (1 + Gg N )) + 0(g 2 N ), (3.25) 



N 

from which we obtain as in (ll.lOp 

/3(g) = -2g - I2g 2 . (3.26) 

Therefore we again find a nontrivial fixed point g* = —1/6 and a critical exponent there 
7i = 2/ f3'(g*) = 1. Although jx is worse than that obtained in [5], they are still compa- 
rable with the exact results g* = —1/12 and 71 = 5/2. There may be several reasons why 
our RG becomes worse than in [5] for the two-dimensional gravity. First we notice that 
our formalism respects the rotational symmetry SO (3). It is evident that this property 
fits in field theory on the (fuzzy) sphere, but not in the two-dimensional gravity describing 
random surfaces. Another reason could be that in our formalism we have dropped con- 
tributions from the antipode matrix, while in [5] all O(g) contributions in the systematic 
1/N expansion are taken into account. 



4 Conclusions and discussions 

In this paper we propose a new large-iV RG based on the fuzzy sphere. It has a nice 
analogy with the usual RG in some aspects such as the locality of the RG transformation 
and the derivative expansion. It also reveals the interesting features in matrix models as 
the appearance of multi-trace operators, the antipode transformation, and renormalization 
of noncommutativity. Fixed point analysis gives consistent results for the Gaussian fixed 
point of field theory on the fuzzy sphere, or the NCFT on the two-dimensional plane. A 
nontrivial fixed point is also found and their properties are discussed. There is a critical 
line if the radius of the sphere is large. Our RG equation provides an expression for the 
critical line which well agrees with the numerical simulation. A comparable value of the 
critical exponent in the two-dimensional gravity can be also obtained. 

Since it seems straightforward to apply our formalism to multi-matrix models, it is 
desirable to apply it to the large- N limit of more interesting models, related to string 
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theory, like the Chern-Simons type matrix models as in e.g. [12j l2Ql IZE]. Our formalism 
is expected to reveal some aspects of them, in particular, on their universal properties. It 
is also anticipated that our large- iV RG also sheds light on the large- N limit of supersym- 
metric gauge theories on higher dimensional sphere that have attracted much attention 
recently. It would be likely that in order to clarify vast universality class of large- N gauge 
theories suggested by the large- N reduced models [22], we would need formalism based 
on the RG like ours. We mention that it would be also interesting to examine a relation 
between our large- iV RG and a master field approach on the fuzzy sphere [23] . 

Before closing the paper, we discuss two issues which would be considerable for further 
understanding of the matrix RG method in the present paper. 

Fuzzy torus 

One may think that our large- N RG is also available on the fuzzy torus, where the algebra 
of functions is more tractable. In fact, a.n N x N matrix <f>N can be expanded analogously 
to the Fourier expansion as (assuming odd N) 



where U and V are the standard 't Hooft matrices with rank N satisfying UV = oo X VU 
with u = exp(2iri/N). Then each vertex of a field theory on the fuzzy torus has a 

factor imposing the momentum conservation modulo iV: S^ N \^ r n ( {^)5^ N \Ylr n 2' 1 )- m 
our large- N RG, it is natural to integrate first over (4N — 4) modes with the maximal 
momentum of n% or n,2, namely, over ni „ 2 with |m| or |rz-2 1 — {N — l)/2, and then 
to rewrite the resulting action as an (N — 2) x (TV — 2) matrix model. However, in 
the (N — 2) x (JV — 2) matrix model the momentum conservation should hold in each 
vertex modulo N — 2. For example, if a vertex consists only of in-modes with |ni|, |n 2 | < 
(N — l)/2, it remains intact in the RG, and hence there the momentum is conserved 
modulo N, not N — 2. This serious problem would be bypassed by considering the large- 
iV RG not by N — >■ iV — 2, but by iV — » N/2 as the block spin transformation. We will 
report our study in this direction elsewhere. 

Antipode, UV/IR mixing, and noncommutative anomaly 

Another interesting aspect in our formalism is the appearance of the antipode transfor- 
mation. Here we point out an interesting connection between the antipode matrix and the 
nonplanar diagram. If we draw nonplanar diagrams by the standard double line notation, 
we find that they resemble annulus diagrams appearing at one-loop of open string theory. 
They have two edges from which external lines emanate. Our intriguing observation up to 
0(g 2 N ) is that if we identify matrices brought by external lines from one edge as "usual" 
ones, external lines from another edge provide the antipode matrices. Namely, there is a 
connection between the Feynman diagram and the positions on the (fuzzy) sphere. We 



jV-l 




(4.1) 
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hope that further analysis of nonplanar diagrams makes clear the origin of the antipode 
transformation. 

Related to this, in [T4] it is claimed that such antipode degrees of freedom trigger the 
UV/IR mixing |17j . After [13] , it is rigorously argued in [19] that there is no UV/IR 
mixing on the fuzzy sphere for finite N based on the one-loop effective action including 
a nonplanar diagram, and the IR singularity seen in field theory on the noncommutative 
plane is reproduced as the large- N limit of a term arising from a nonplanar diagram, that 
the authors call noncommutative anomaly (it is regular for finite N). As seen below in 
detail, the appearance of the antipode matrix reported in this paper can be understood 
as a part of such a noncommutative anomaly for finite N. 

Here we make a detailed comparison of (V^ fP ^ between that in this paper and in [T9] . 
The authors of [19] use the following formula in calculating (V^ p ^ 

Pl U_Jl\. (4 . 2) 



LLJj 2L l \ 2L\ 

The noncommutative anomaly for finite N is obtained as the difference between (V^ fP ^ 
and \V P \ with the internal mode integrated over the whole region < J < 2L by use 
of this formula. On the other hand, our antipode matrix originates in (V^ p ^ with only 
the J = 2L mode integrated and with (12.671) applied. The term with J = 2L in (14.21) 
indeed reproduces the first term0 in the formula (12.671) . Therefore, the antipode effect of 
(V£ iP ) Q in this paper can be regarded as the maximum momentum part of the noncom- 
mutative anomaly discussed in [19] . It is then preferable that the other appearances of the 
antipode matrix can also be interpreted in a similar way. If we can repeatedly integrate 
out the highest modes, the total antipode effect will eventually, in principle, reproduce 
the noncommutative anomaly for finite N. Our RG approach would be helpful toward 
such a problem, as discussed next. 

Since as shown in [19] the noncommutative anomaly develops the IR singularity in the 
large- N NCFT limit (I2.52p and triggers the UV/IR mixing, the above observation implies 
that around a fixed point corresponding to the large-iV limit in our RG, interactions 
involving the antipode matrix would be relevant and that analysis of them there would 
reveal universal nature of the noncommutative anomaly and the UV/IR mixing. More 
generally, an important question from our study is on existence of nontrivial theory with 
the antipode matrix. According to the spirit of the RG, we are led to consider a new 
matrix model like 

Sn = Ntr N (—[Li, 0] 2 - a[L h 0] [L u A ] + + p\m 2 ^ A + ■ ■ -V (4.3) 

and look for a nontrivial fixed point of our large-iV RG. For example, the fixed point 
analysis of m 2 is expected to be useful for understanding the UV/IR mixing. 



16 Note that (I4.2[) and (|2.67p are different approximation formulas even when J = 2L, and the "derivative 
corrections" in (|2.67l) are not visible by use of (|4.2p . 
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Finally, we comment on the observations made in [T7J[24]. In these papers the UV/IR 
mixing is understood as the exchanges of light degrees of freedom which are called closed 
string modes. This nicely fits with stringy interpretation of a nonplanar diagram. Recall 
that the 1-loop correction of the scalar field 2-point function is UV divergent in the 2- 
dimensional continuum theory and that if we try to regularize it via the noncommutativity 
as in the fuzzy sphere or in the noncommutative plane, a nonplanar diagram in general 
yields the noncommutative anomaly, or the UV/IR mixing, in the large- N limit. Then 
the IR singularity there can be regarded as the effect of the propagation of the closed 
sting mode. Therefore it is natural to expect that even for finite N a nonplanar diagram 
of the 2-point function in the matrix model on fuzzy sphere has a counterpart of such a 
closed string mode. In fact, the formula (12 .68j) can be written as the exchange of particles 
with propagator 

following the expression in [24] . Here identification is made as 

l(l + l)=p 2 p 2 , <^p 2 =po P , (4.5) 

and p is the radius of sphere. The spectral function X(p 2 ) in (24] has a uniform distri- 
bution which vanishes below l/(pA) 2 , where A is a UV cutoff of Wilsonian RG in NCFT 
and p is replaced by a parameter with dimensions of squared length in the reference. On 
the other hand, the spectral function in (14 .4p is the delta function. This is because we 
calculated a one-step RG flow N — > N — 1. It is shown that matrix models on the fuzzy 
sphere are free from the IR divergence for finite A" [19]. In fact, this feature can be also 
seen in ( |4.4p . The spectral function does not contain massless modes unless A" is infinity. 
An IR regularization of the closed string picture is realized in this way. 
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A Useful formulas of 3?r/-symbols 

In this appendix we enumerate useful formulas for our study from [15] . The Latin and 
corresponding Greek indices denote angular momenta and magnetic quantum numbers 
(projection of angular momentum), respectively. For example, the magnetic quantum 



30 



number a runs —a < a < a for the angular momentum a. In the following we define 

, , , _ J 1 \a — b\ < c < a + b , , 

^ ^ ^ otherwise \ ■ J 

namely, {abc} does not vanish only if a, b, and c satisfy the triangular conditions which 
are symmetric under the interchange of variables. 

A.l Identities 
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A. 2 Asymptotic formulas of 6j-symbol 

If a, b, c 3> / and / is an arbitrary integer, 

a b c 1 (-iy+b+c+f 



•) ( _l)a+b+c+f 

i f > - , =P f (cos 9), (A.9) 

bafj v /(2a+l)(26+l) A 



with 



cos 9 = ■»(■»+ D + K*+l)-^+l) (A . 10) 

2^0(0+1)6(6+1) 



If i? ^> 1 and a, 6, c are arbitrary, 



\c+d+e 



\d + i2e + i2/ + i2/ - ^2i2(2c + 1) < f ~ e) ^ ' 1 J 

where C c ^ ahj3 is the Clebsch-Gordan coefficient^^ 

To investigate the asymptotic behavior of the 6j'-symbols, it is sometimes useful to 
use the explicit expression due to Racah, 

{ a f c } = g ' C)A(A ' 6 ' C)A(Q - B ' c)A(a ' 6 ' C) E ( ~ 1) /( t t) + » ( A - 12 ) 

/(t) =(t - m)\(t - n 2 )\{t - n 3 )\(t - n 4 )!(mi - t)\(m 2 - f)!(m 3 - t)\ , 

a 1 \ _ (x + y - z)\(y + z - x)\(z + x-y)\ 
i\{x, y,z) — . . , 

(2 + y + z + 1)! 

ni =A + 73 + C , n 2 = A + b + c, n 3 = a + B + c , = a + b + C , 

m 1= A + a + B + b, m 2 = B + b + C + c, m 3 = C + c + A + a, 

where the sum of t runs in the region where no argument of the factorials takes negative 
values, namely, 

max (rii, TI2, n 3 , n 4 ) <t< min (m l5 m 2 ,m 3 ) • 

Several asymptotic formulas are obtained by evaluating the factorials by the use of the 
Stirling's formula. 

B Useful formulas of the fuzzy spherical harmonics 

In this appendix useful formulas of the N x N matrix T/ m defined in f)2.12p are presented. 
We set L = (N — l)/2 as in the main text. (1A.4I) . ( 1A.5I) and ( 1A.6I) give formulas of the 

trace of T; m 's 



3 / 

trjy (T hmi T hm2 T hmz ) = Nl Y[(2k + l)l(-l)^+Ef=i^ I 

i=i ^ 



h h h \ \ h h h 
mi 777,2 J I L L L 



17 



We have supplied the formula presented in [15] with a necessary phase which was also added in [19] 
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trjv (Ti 



N 2 X\{2U + lt^- l Y m ( 2l + l ){ 

i=l Ira ^ 



mi m I \ —m m2 




N 2 Y](2k + i) l 2(-iy^y(-iy~ m (2i + i)( h h l ) ( 1 h 1 

i=l lm v / \ - 



Notice that (IB .2j) can be also derived easily from (IB.ip and ( I2.14p . and that U (i = 
1 ~ 4) on the right-hand side of flB.2|) can be cyclically interchanged, which reflects the 
cyclic symmetry of the trace. From (12 . 1 3[) and (IB.lj) we can read off a fusion coefficient 
F hmi i 2 m 2 hm3 defined by 

Tl im \Ti 2m2 ^ Fl\rn\ l 2 m 2 -^3"i3' (B'^) 

as 

F hmikm >™ = TV* n (2ii + l)J(-l) 21 +£!.i '<+">" £ jy { ^ } . (B.5) 

As an application, by using the asymptotic form of the 6j-symbol for L ^> 1 in ( 1A.11I) 

olcl (-l) c 
£££J ~ v/2L(2c+f) 

and 



1) 2L ,,15, (B.6) 



( h h h \ _ / i\< 3 +w 3 +2ii -*- fyl 3 m 3 (B7) 

(IB .41) and ( IB. 51) reproduces the fusion of the usual spherical harmonics [15] 



v v - V * / ( 2/ i + 1 )( 2 ^2 + 1) ^30 w 3 m 3 v m 8^ 

-f limi -f l 2 m 2 — Y 2/ + 1 hOhO hm 1 l 2 m2 1 hra- i - K^- ) 

Hence the structure constant of the fuzzy spherical harmonics 

Plimi > -^2mJ — fhmx hm? 23m3-^l 3 m 3 > 

^3^3 

f / t \ m,3 / en h — ms _ 771 i 3 — m 3 \ m q\ 

Jhmi l 2 m 2 £3)713 V J-J y-^iimi Z 2 «i2 r l2m2limi ^ 5 V- 0,3 / 

also coincides with that of the usual spherical harmonics in the large- iV limit. 
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Furthermore, by using the fusion f ]B.4j) and flB.5j) repeatedly, we find 
T 2Lm T hmi T 2L . m =N(AL + ( 2l ' + 1)7(2/1 + l)(2/"+l)(-l) h+r+m ' +m " 



I'm' 
l"m" 



X 



h V 2L \ ( V 2L I" \ ( 2L h V \ ( V 2L I" 



mi —rn m J \m —m —m 



L L L L L L 



Tl" m " ■ 

(B.10) 



Taking the sum J2 m =-2L 

2L 



l) m of both sides of this equation and applying (1A.2j) . we have 

( — ^) m T2LmTl 1 m 1 T2L-m 

l) 2L+h N(2N - l)T hmi ^(-l) r (2/' + 1) [ 

i' 



m=-2L 



2L h I' 
L L L 



B.ll) 



where the summation over I' can be also taken by utilizing (IA. T[) to yield 
(-±) m T2LmT hmi T 2L _ m = N(2N - 1) | L L L L ^ ] (-l)' 1 ^. 



(B.12) 



m=-2L 



This equation tells us an interesting fact that making the similarity transformation with 
respect to T 2 L m and taking the summation over m is essentially equivalent to perform- 
ing the "antipode" transformation Ti m — > (—l) l Ti m , namely transporting a field to the 
antipode on the sphere. As a corollary, we find that 



2L 



( — 1 ) m T 2 L m T 2 L -m = {2N ~ 1)1 



N- 



(B.13) 



m=-2L 

where we have used T 00 = In and 



L L 
L L2L 



1 

N' 



(B.14) 



(IB. 121) and (1B.13j) are of great help in calculating our RG. 



C Calculations of 0(g 2 N ) 

In this appendix we evaluate contributions of 0(g N ) in the large- RG shown in the 
second line in (I2.4ip . First we consider (V 2 2 ) c that contributes to vertex correction and 
then turn to the other terms giving rise to mass correction. 
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C.l Vertex correction 

V 2 = V 2 P +V 2 NP as in (I2.36P and the only possible planar diagram in (V 2 2 ) c arises in (V 2 p2 

calculated in ( I2.74p . Furthermore, in the low energy regime where the momentum / of 
(p 1 ^ is much smaller than the cutoff 2L, we can perform "derivative" expansion of it as in 

(12.761) to get (12.781) . Similarly we can evaluate the nonplanar diagram in (v 2 p2 ^ given 

as the second term in (12.701) and obtain a similar result to (12.741) 

2 



Vn 



P2 



nonplanar 



I in 



2L 2L I 
L L L 



-i)W 



1)111 Ti m ) trjv 



Here we note that 



:-i)W 



■t 



I in 



tr 



N 



A 2 



T 



l in 



Im 

(C.l) 
(C.2) 



which can be readily shown directly from the expression of the trace of three T; m 's given in 
( IB . 1 P by using the well-known property of the 3j-symbol, or from the proposition proved 
in the appendix [D] Then in the low energy regime the same derivation as in (j2.75p ~( l2.78p 
leads 



V, 



P2 



nonplanar 



P 2 N {2N - l)N 



i2 mtij4^\ ^ 

' 2N' 



Al\ I 1 



N 2 



(C.3) 



On the other hand, from ( 12.361) we obtain 

{vfv?-) c 

-\ E 



mi,rn2,m' 1 ,m2 



lout i out 



lout J,out 



'2L n, 1 ' ''21. n,,' 2 / \4 > 2Lm 2 ( f > 2Lm' l / + yf > 2Lmi ( l ) 2Lm' 1 j \V2Lm 2 ^2Lm,' 2 f ( . 



x trjv (0 in T 2Lmi T 2Lm2 j tr N (0 in T 2Wi in T 2 w 2 ) 
P 2 N (-l) m+m 'tr N (0 in V 2im T 2fcm ,) tr N (0 in T 2L _ m ,0 in T 2L _ m ) , (C.4) 



where the two kinds of the contractions yielding nonplanar diagrams turn out to be 
equivalent from the cyclicity of the trace. By expanding in2 as m2 = Yli m ^hn Ti m and 
using the trace formula in (IB.ip and in (1B.2j) . we have 



P 2 N {2N-lfNl-£^) m+m ' E C2< 



m,m' 



(( r ,m r )6A; n 



35 



x (-1 



\2L+h 



1 + 1 

r=l 

^(-ir"(2r + i) 



i Mi 2L 2L 



l"m" 



mi in m 

l 2 2L I" 
m 2 —m —m! 



h 2L 2L 
L L L 

I" h 2L 
m" m 3 —m 1 



l 2 2L I" 
L L L 



I" l 3 2L 
L L L 



(C.5) 



Here we can take the summations over m, m', and m" according to (IA.4[) . while the sum 
on I" can be also taken by (1A.8|) . We thus get 



E 

(i r ,m r )6A; E 
r=l~3 



,in(2) ,in ,in 



\[{2k + l)H-l) 1 



i=l 



h h h 
mi ni2 "^3 



h 2L 2L 
L L L 




(C.6) 



We next use the cyclic symmetry of the 9?-symbol and utilize flA.8j) once again. Then the 
9j-symbol can be rewritten as 







'h hh) 


\LL2L 


H 


2L L L } 


[l L 2L j 




2L L L) 



h h h 
L L x 



2L L L 

l 2 x L 



2L L L 
x h 2L 



(C.7) 

Note that in contrast to I" in (IC.5j) . here the selection rules of the 6j-symbols impose 
L < x < L + min(/i, / 2 ) and x G Z. Thus setting x = L + m, we have 



min(Zi,Z 2 ) 



-1) 2L J2 (N + 2m) 



m=0 



h h h 
L L L + m 



2L 



I2 L + m L ) \L + m li 2L 



2L L L 



(C.8) 



Three 67-symbols above can be explicitly evaluated via the Racah formula (IA.12P and 
it can be further approximated by the Stirling formula in the low energy regime l r L 
(r = 1 ~ 3). For instance, 

2L L L 
I L + m L 



(l + m)\ / ((2L)!) 2 (2L + m)!(2L -m)!(4L + m + 1)! 



<m! y (I - m)\ y (4L + 1)!(2L - l)\(2L + I + 1)!(2L + m - l)\(2L + m + I + 1)! 
(-l) m l(l + m)\ 1 



(/-m)!2Lf +1 



1 



8/(7 + 1) - m(m - 5) + 



16L 







1 

1? 



(C.9) 



where we note that terms with higher m are suppressed as L 1 2 when L ^> 1. We can 
evaluate the last 67-symbol in (1C.8I) in a similar manner, 



2L L L 
L + m I 2L 



(l + m)\ Vg 
/^T y (i-m)! (4L)f +1 
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2/(/ + l) + m(5m + 7) + 6 
16L 



and this is also suppressed as L -1- ^ '. Therefore, in the 1/L expansion, it is sufficient to 
consider m = and 1 in the summation. As for the remaining 6j-symbol, we can relate 
m = 1 one with m = case for I, < L (i = 1 ~ 3) 



h h h 
L L L + 1 



h{h + 1) - h(k + 1) - h(k + 1) / h h h 

L L L 



2^h(h + l)l 2 (l 2 + l) 



(CIO) 



by use of the asymptotic formula (IA.110 and a recursion relation for the Clebsch-Gordan 
coefficient 



°al6-l 



c(c+ 1) - a(a + 1) - 6(6+ 1) 



aOMT 



2^0(0 + 1)6(6+1) 

Combining these, we obtain 

min(£i,£ 2 ) 

(AT J- 9ml J 

L L L + m\ \ l 2 L + m L ] \L + m l\ 2L 
1 



(C.11) 



dn(i 1 ,i 2 ) , 

^ (N + 2m) J 

m=0 ^ 



h h h 



2L L L 



2L L L 



ir h h h 



y/2N \L L L 



1 - 



N V 4 



--Z^ + 1) + -Z 2 (Z 2 + 1) + -z 3 (z 3 + 1) - 7 1 + o 



1 

(C.12) 



Plugging ( jC 121) into (IC.6I) . together with (12 .75 p . we see that the remaining factors are 
exactly of the form of ( IB. II) and finally obtain 



= P%(2N-1)N 



tr 



iV 



2A 



"2iV triV 







AT2 



(C.13) 



The calculation of (v 2 Np2 ^ is similar. The two ways of contractions again give the same 
result as 



p2 



1 T2Lm3<p m T2Lm 4 ) ^ T N 



—m,4 



p2 
N , 



. 2 EM)- 

m3,m4 

2 

(2A - 1) 2 A 4 II ((«r + VHtmM + l) 1 ^^) (-I] 

l\ l 2 I \ ( Zi Zn I 



h+l' 2 



{l r ,m r ),(l' r ,m' r )eA in r=l 
r=l,2 



y(-iym (2Z + 1} ( «i fa z W z; z 2 z \ 

^— ' \mi m 2 —my \m 1 m 2 ml 



r Zi z 2 z > 






'.LL2L 


w 


L L 2L J 


[l L 2L j 




L L 2L) 



(C.14) 
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where we have used the formula for the trace of four T/ m 's in terms of 9j-symbols flB.3j) 
instead of 6j'-symbols (IB. 21) . Then we apply a similar evaluation from flO.Zj) to (IC.12P 
to two 9j-symbols in this expression. A lengthy calculation as above results in the last 
equation in (12.79p . 



C.2 Mass correction 

(V 1 2 ) c would give rise to a 6 vertex. But we shall see that it is negligible in the low energy 
regime. In fact, from ( I2.36p . (V 1 2 ) c becomes 

( V i) c = E WmOo^ (<P in3 T 2Lm ) tr N (f^T 2Lm ) 

m,m' 

=P N ^(-l) m t rjv (0 in3 T 2Lm ) tr N (0 in3 T 2L _ m ) . (C.15) 

m 

■ 3 

However, when the momentum I of m is much smaller than 2L, that of m cannot be 
equal to 2L and the traces above vanish by (I2.13p . Hence in the low energy regime of 
interest in the RG we can neglect this contribution. 

The other remaining terms of 0(g 2 N ) will contribute to quadratic terms of m . First 
let us consider (V 2 ) c . Using the trace formula (IB.2p . we have after contractions 



{V 2 ) c =P" N (2N - l) 3 iV 4 Krn^VW+m^+l) 

(.h,mi),(l' 1 ,m' 1 )eA in 

h 2L l\ ( l[ 2Ll'\( I 2L2L\ ( V 2L 2L 
L L L L L L L L L \\LL L 



£(2Z + 1)(21' + 1){ 
i,i' ^ 



X 

mi ,~,m3 m,m 

l\ 2L I' \ ( V 2L 2L 
x . i ( 



l x 2L I \ f I 2L 2L 

mi rriA m J \ — m m 2 



m'i — m 2 m! ) \ —ml —m^ —m^ 



+ (permutations of m 2 ~ m^j 



(C.16) 

When lx, l[ <C L, the first two 6j-symbols impose / = 2L — m, I' = 2L — n with m <l\, 
n < l[. Then the Racah formula (IA. 12[) and the Stirling's formula show that the large- L 
behavior of the third 6j-symbol as 



2L-mQ±(r,\± _ / Q \ 3L 



2L-m2L2L\ (-l) 2L " m 34 (2tt) i 



3 



, , , ^- L^-i ^ , (C.17) 

namely, it is exponentially suppressed for L> 1. Since it is easy to see that other factors 
in the above equation are bounded at least by polynomials of L, we conclude that (V 1 2 ) c 
only gives exponentially small contribution in the low energy regime. Furthermore, it is 
easy to show that (ViV3) c = by the momentum conservation. Thus only (V^V^ provides 
a nonzero contribution even at low energy. From (I2.36P we first find that 



(V 2 p V 4 ) c = P 3 N £ (-l) mi+m2+m3 tr w (0 in2 T 2Lmi T 2im2 ) 



mi,- ,ni3 



38 



x (^> T N (T2L-m 2 T2Lm 1 T2LmT2L-m) + ^ T N (?2L -mi^L m-^%L m^2L -m) 
+trjV (T2L-m 2 T2LmT2Lm 1 T2L-m) ), (C.18) 



where the first two terms correspond to planar diagrams, while the last term to a nonplanar 
one. The formulas (1B.12j) and (IB.13j) enable us to rewrite this as 



(V 2 p V 4 ) c = P" N (2N - l) 2 N ^2 + {-lf L I 



LL1L 
LL1L 



tr N U in ' . (C.19) 



In the second term that comes from the nonplanar diagram, the 6j-symbol is again expo- 
nentially small as 

LL2L) /27T„_ ]L _ 2 



L L 2L $ ~ )/ L ' (C20) 

for L>1. Therefore for L ^> 1 the planar diagram gives 

(V 2 P V A ) c ~ 2(2iV - l) 2 iVtr w (> 2 ) . (C.21) 

Finally, we evaluate (V^ P V4) . This can be carried out by using the technique ex- 
plained so far. We quote only the result, 

{V 2 NP V 4 ) c =N(2N - 1) 2 P^ tTN _ ^ m^j j + . . /j (C 22) 

Therefore, this term involves the antipode fields. 

D Antipode transformation and ordering reverse 

In this appendix we prove a proposition that reveals an interesting connection between 
the ordering of matrices inside the trace and the antipode transformation. 
Proposition: 

Proof : 

The n = 1 case is trivial. The n = 2 and n = 3 cases are also obvious because of the 
orthogonality ( 12 . 1 3[) and the explicit expression of the trace of three generators ( IB. 1|) with 



the symmetry property of the 3j-symbol. Assuming ( ID.lj) for n = 1, • • • , k, from (I2.15P 



f k+l \ /fc-1 \ 

tr -v ( 11^ = n [U^Tim tr N (<t>M +1 Tl 

\i=l / Im \i=l / 
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^5>m (-i)^n^- tr ^ ((- i )' T / m ^+i < 

Im \ i=l / 

tTN ( II fa-iTim ) tr N (cj) k+l (j) k T ] lm 



Im \i=l 

fe — 1 \ /a-;-.. 

tr N [ <p k +i<Pk 

i=l / \i=l 



[J 4>h-i = , r v I I </>k+2-i , (D.2) 



thus (ID.lj) holds for n — k + 1. This completes the induction. 
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